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Abstract 

By using Laplace transformation we developed an approximate so- 
lution to describe neutrino oscillation probabilities in arbitrary density 
matter. We show that this approximation solution is valid when matter 
potential V satisfy V < Am'^/2E and / VL < 1, where L is the length 
of the neutrino oscillation . Thus, the formula is useful for propagation 
of the solar or supernova neutrinos with terrestrial matter effect. 
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1 Introduction 



By now, we have solar, supernova, atmospheric, reactor and accelerator 
neutrino experiments to determine the mass splittings and flavor mixings 
in neutrinos. In order to calculate the survival and conversion probabilities 
of neutrinos passing through the Earth, the matter effect ^ must be taken 
into account. Previous work on this subject inculdes exact and approximate 
expressions in the case of constant matter density [2]-|17j. linear density |S1 
and exponential density [7j. As for the case of arbitrary matter density, 
approximate solutions have also be presented by Akhmedov ^S] , Peres 
and loannisian j2Uj. j21j. 

Motivated by loannisian's paper j^Uj, we derived an approximate solution 
to MSW equation by using Laplace Transformation. The form of the formula 
in this note is so simple and transparent that it is no difficult to calculate 
up to the high order terms. So, with this formula, we could simplify the 
numerical calculation considerably. We also show that the solution is valid 
for the case of low energy neutrino such as solar and supernova neutrinos, 
and the approximation is effective when the baseline length is not too long. 

The outline of the paper is as follows. First, we investigate the general 
approximate solution for electron neutrino survival probability with 2 fla- 
vors. Then we compare our formulas with exact numerical results in some 
special case: uniform matter density and linear matter density. And finally 
we discuss some related issues. 

2 General approximate formula by using Laplace 
Transformation 

We consider the case of two-flavor (electron flavor and the effective flavor 
ly^ - a linear combination of and i^r) neutrino oscillation for simplicity. 
The neutrino mixing between flavor-eigenstates and mass-eigenstats could 
be expressed as 



and mass eigenstates respectively. And the lepton mixing matrix is given 





with Uf = (fe(x), Uxix))"^ and 



(z^i(x), iy2{x))'^ are the flavor eigenstates 



by 
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In order to find the neutrino oscillation probabilities in matter, we have to 
solve the Schrodinger equation for flavor eigenstates 



; — 

dx 



U,j:{x) 



H{x) 



Ue{x) 



VAX) 



And the effective Hamiltonian for propagation of neutrinos in matter 



U 



nil 



A{x) 







(3) 



(4) 



where A{x) = 2^/2GpNf,{x)E is the effective potential term. With Gp is 
the Fermi constant, Ne{x) stands for the electron density in matter at point 
X, and E is the neutrino beam energy. We take 



3.8 X IQ-'^eV^ 



YeP 



1 GeVj 



(5) 



^2.5g/cm^ 

Notice that the sign of the matter potential is positive for neutrinos and 
negative for anti- neutrinos. 

Applying Laplace transformation to the Schrodinger equation, we get 



IS 



+ 



C[Mx)] J [ V,{0) J 

Ai cos^ + A2 sin^ d 
- Ai sin 9 cos 6 + A2 sin 6 cos 9 

C[V{x)Ue{x)] \ 





- Ai sin 9 cos 9 + A2 sin 9 cos 9 
Ai sin^ 9 + A2 cos^ 9 




Here we define Ai = mj/lE, A2 = mi/2E and V{x) = A{x)/2E. To 
find the relation between I'eix) and V{x), we rewrite this equation as 



CWeix)] 

C[i^x{x)] 



is — Ai cos^ 9 — A2 sin^ 9 Ai sin 9 cos 9 — A2 sin 9 cos ( 
Ai sin 9 cos 9 — A2 sin 9 cos 9 is — Ai sin^ 9 — A2 cos^ 9 



MO) 

MO) 



+ 



c[v{x)Mx)] 




It is straightforward to set the initial condition as z^e(O) = 1 and i^xi^) = 
when we calculate the electron neutrino survival probability Pee and con- 
version probability Pex- Thus Eq. © could be expressed explicitly as 



c[Mx)] 



COS 



(s + iAi 



+ 



sm 



> + iA2) 



(6) 



-1 



(7) 
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cos^ 6 sin^ ( 



and 

1 1 



C[V{x)i,,{x)l (8) 



s + i/S.1 s + iAi, 
-z sine cos (—1- C[V{x)ve{x)]. (9) 

Then we apply Inverse Laplace Transformation to Eq. (jH)) and @, and 
we arrive at two integral equations of the oscillation amplitude of Ve and Vx'- 

ue{x) = cos^ 06-^^1^ + sin^ 0e-^^2^ 

-i£ dy (cos^ 0e-^^i(^-^) +sin2 ^g-^^^^^"^)) y(y)i/e(2/)(10) 



= sin0cose(e-^^i^ 

-i sin cos £ dy (g-^^^^^"^) - g-^^^Cx-y)^ V{y)iye{y). (11) 

Notice that the third term in the right side of Eq. H10|) and Ijlip is a con- 
volution integral. We will deal with Eq. (|lfl|) first, and it is clear that Eq. 
(jllf) could be solved straightforwardly with the result of Eq. (fTU)) . 
In order to get the approximate solution to this equation, it is convenient 
to define an operator 

K{ue{x)) = -i £ dy (cos^ 0e-^^^(^-S') + sin^ 0e-i^2ix-y)^ V{y)iye{y) 

= -icos^ee-'^''' r dye'^^yV{y)v^{y) 
Jo 

-iWOe-'^^y r dye'^^yV{y)ve{y)- (12) 
Jo 

Thus, Eq. ()12p could be rewrited as 

ve{x) = cos^ 0e-*^i^' + sin^ ^e'^^^x ^ K{ve{x)). (13) 

Thus, if II K{ue{x)) II < 1, the solution of Veix) could be expressed in a series 
expansion form: 

veix) = (l + K + K^ + K'^ + -->j (cos^ 0e-^^i^ + sin^ ^e'^^^^) . (14) 
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Inserting Eq. H13|) into Eq. the oscillation amplitude of Vx 

ux{x) = sin^cose(e-^^i^ 

-i sin e cos dy (g-^^^^^-?^) - e-*^2(^-?^)) F (y)zye(y). (15) 

A straightforward calculation leads to the Ue survival probability P^^^^e = 
|t'e(a;)P and conversion probability Py^^y^ = \vx{x)\^- We note that other 
oscillation probabilities such as P^^-^u^ or P^^^y^ could be obtained with 
the same method but different initial condition {[Vf>,Vx^ = [0, 1]"^)- 



3 Numerical discussion and testing of accuracy 

Eq. (|14jl in the last section is a general formula. In this section, we discuss 
the qualitative behavior of this formula. First let us consider the case of 
constant matter density {y{x) = V). And with the result we can derive a 
raw applicability condition for this formula. 

Eq. (|14|) could be expanded in orders explicitly as the matter potential 
is a constant V{x) = V . The first three orders of the amplitude of electron 
neutrino oscillation Ve{x) are 

= cos^ 0e-*^i^ + sin^ 06"'^^^ (16) 



z,g(x)W = -zcos^eyxe-*^i^-isin^0yxe^*^2x 

91/ 

-sin^^cos^e— ^— — (e-^^2x _ -iAix) ^^^^ 



ve{xf^ = -cos^e^e-^^i^-sinS^^e-*^^" 

-3i sin2 e cos" e-'^^" + 3i sin^ 9 cos^ 6 ^^^""^ e"^^^" 

Ai - A2 Ai - A2 

+-sin2 2ecos20-— — (e-^^2x _g-*AixN /^^g) 
4 (Ai-A2)2^ 

We can see that there are two expansion parameters in Eq. ()16() - (|18p : Vx 
,which has appeared in ^H] and ^^^^^ . If the number density of the 
electrons = 2.26 and x = DEarthi the diameter of the earth, we 
could find that Vx ~ 0.4. Since in the case of earth-induced matter effect, 
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X < Dearth, we could estimate that Vx < 0.4 thus it could be proved as 
an effective expansion parameter. Another expansion parameter ^^^^^ is 
similar to the e(x) in loannisian's paper and which would be small 
enough ( ^^^^^ < 0.6 if < 50MeV |U)for approximation in solar and 
supernova neutrinos (low energy). 

The oscillation probability up to the zero-th order (we set Ai = 0): 

PHe = |-e(x)Mp 

= 1 -sin2 26'sin2(^x), (19) 

first order: 

PHe = ke(x)m + ^e(x)Wp 

= 1 - sin^26lsin^(^a;) 

Vx V 
— — sin^ 29 cos 29 sin A2X + sin'^ 29 cos 20(1 — cos A2X) , 



(20) 



second order: 



pi% = |^,(^)[0]+^^(^)[l]+^^(^)[2] |2 



1 -sin^ 26'sin2(— x) 
^ 2 ^ 

Vx V 
— — sin^ 29 cos 29 sin A2X + — sin^ 29 cos 26'(1 - cos A2X) 
2 A2 

sin^ 49 cos Aox 

16 

^ ^- sin^ 9 cos^ 6* (4 sin"^ 6^-12 sin^ 9 cos^ 9 + 4 cos'' 6*) sin Agx 



A2 

(21) 



+ ^ sin^ 9 cos^ 6* (6 sin^ 6* - 20 sin^ 9 cos^ 6^ + 6 cos^ 9) (cos A2X - 1 



Inserting Eq. (|T6|) - (fT7j) to Eq. (fTT|) . we get the amplitude of Uxix) up to 
the first order 



z/^(x)[o] = sin0cos0 (e-*^i^ - e" 



(22) 
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and 



-i sin e cos^ eVxe''^^'' - i sin^ 6 cos 6iyxe-^^2X) 

+- sin 20 COS 261 . (23) 

2 Ai — A2 

And from Eq. (|23p . we get the oscillation probability (omit the second 
order term) 



= sin^20sin^( — -x) 



Vx V 

+ sin^ 26 cos 20 sin A2X — — sin 20 cos 20(1 — cos A2X). 

2 A2 

(24) 

We find pjlle + Pl'ix = 1, which could serve as a cross check of the 
formula. 

Since in the case of constant matter density, the convergency condition 
of this approximation is Vx < 1 and ^^^^^ < 1, we can arrive at a raw 
convergency condition for arbitrary density: 

VmaxX < 1, (25) 

and 



Ai - As 



< 1. (26) 



It is transparent that both of the two conditions could be widely satisfied if 
we investigate the terrestrial matter effect of low energy neutrino. 

In order to examine the reliability of this formula, we compare the elec- 
tron neutrino survival probability Pe-»e obtained from the zeroth order cor- 
rection, first order correction and second order correction with the exact 
numerical result. In Fig. and we observe that the solution to the 
electron neutrino survival probability Pe^e is a good approximation already 
when the length of neutrino propagation L < 6000 km and neutrino energy 
E < 20MeV. As for the case of L > 6000 km and E > 20MeV, we have 
to calculate the approximate probability up to the second or third order. 
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When the potential is a linear function V{x) = a + bx, the formula is 
also useful. According to Eq. (|14() we can arrive at an approximate solution 
up to the first order. 

i/e(x)[°l = cos^ ^e-^^i^ + sin^ 6'e-^^2^ (27) 

Ueixp = -icos'^eaxe-'^^'' -isin^eaxe-'^^'' 

- sin2 6 cos2 0—^^ - e-'^^^l 

-i cos^ - i sin^ 0— e'^^^^ 

2 2 

(28) 

So, the oscillation probability (we set Ai = 0) 

P^e = ke(x)[°l+^e(x)Wp 
9 9 A2 

~ 1 — sin 26* sin ( — x) 
^ 2 ^ 

- ^ sin^ 261 cos 29 sin A2X + sin^ 20 cos 26* ( 1 - cos A2X) 
2 A2 

— ^sin2 2(9cos26'sinA2X. (29) 

Here we calculate a special case: when the neutrino flux pass through the 
center of the Earth, according to 23 , the potential could be expressed as 
four sections of linear functions approximately. 

V{x) = ai + bx < X < xi, 

V{x) = a2 + bx xi < X < X2, 

V{x) = as — bx X2 < X < X3, 

V{x) = — bx Xs < X < X4, (30) 

where ai = 1.14x lO'^^ [eF], 02 = 3.04x IQ-^^ [eF], 03 = 6.84x IQ-^^ [eF], 04 = 
4.94xl0-^3[ey], b = 2.7xlO-'^'^[eV]/[km], xi = 2800km, X2 = 6200/cm, xa = 
9600A;m, X4 = 12000/cm. With Eq. ^ and Eq. (gHl), we obtain 

Pe^e = 1 - sin^ 26 sin2(^x) 

+2(cos^ 6 + sin^ 9 cos A2X4) (/ sin^ 6 sin A2X4 + sin^ 6 cos^ 9) 

A2 

+2 sin^ sin A2X4(— / cos^ — f sin^ 6 cos A2X4 + sin^ 9 cos^ 9), 

A2 

(31) 
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where 

f = aixi + 02X2 - a2Xi + 03X3 - 03X2 + 04X4 - 04X3 + 6x2 (32) 

gi = ai cos A2(x4 — xi) — ai cos A2X4 + 02 cos A2(x4 — X2) 

— a2 cos A2(X4 — Xi) + 03 COS A2(X4 — X3) — 03 COS A2(X4 — X2) 

—04 COS A2(x4 — X3) — ai COS A2X1 — 02 cos A2X2 + 02 cos A2X1 
—03 COS A2X3 + 03 COS A2X2 — a4 COS A2X4 + cos A2X3 + oi + 04, 

(33) 

52 = ai sin A2(x4 — xi) — ai sin A2X4 + 02 sin A2(x4 — X2) 

— a2 sin A2(x4 — xi) + sin A2(x4 — X3) — sin A2(x4 — X2) 
—04 sin A2(x4 — X3) — oi sin A2X1 — 02 sin A2X2 + 02 sin A2X1 
—03 sin A2X3 + 03 sin A2X2 — 04 sin A2X4 + 04 sin A2X3, (34) 

The numerical result is showed in Fig. ((31), compared with the result of uni- 
form density (we set the potential ^ = 3 x 10~^^[eT^], which is the average 
density of the Earth). And we can find that in such case, constant potential 
is not a good approximation. 

Furthermore, the formula of two neutrino species could be generalized 
in a straightforward way to the case of any neutrino species. If there are N 
types of neutrino involved, Eq. Q turns into 

^iMx)] = 7— -ITT + 7— -ITT + ■■■ + 



(s + iAi) {s + iA2) (s + iAN) 

+ , . . . + ••• + 



Xs + iAi) (s + iA2) {s + IAn) 

xC[V{x)iye{x)], (35) 

where 

at = UeiUli, (36) 

and U is the N-flavor mixing matrix. Thus, the solution could be expressed 
as 

ve{x) = (l + K' + K'^ + K"' + -->j (aie-^^i^ + 026"^^^^ + • • • + o^e"^^^^) . (37) 
Here the operator K' is redefined as 

K'{ve{x)) = -i f dy (aie-^^^^^"*') + a2e-^^2(^-^) + • • • + aive-^^^^^^"^)) V{y)ue{y). 

(38) 
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4 Conclusion 



We have derived a series expansion formulation for neutrino oscillation in 
arbitrary density matter. Examining the case of constant density, we found 
that this formula is useful in calculating electron survival and conversion 
probability with a not too long baseline { L < 12000 km ) and low energy 
( £^ < 70 — 80 MeV ). As for some special cases, we can see that this 
expansion up to second order is possible to reach a high accuracy. We 
also show that this formulation could be extended to the case of N-flavor 
neutrinos by adding similar extra terms to the expansion operator K. 
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Figure 1: Electron neutrino survival probability as a function of length of 
baseline L with E = lOMeV and p = Q.Og/cm^ {V ^ 2 x and 
best-fit KamLAND data [HI parameter values: Am^^^ = 7.1 x 10^^ eV'^ and 
&soi = 32.5°. 
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Figure 2: Electron neutrino survival probability as a function of energy E 
with baseline length L = 7400 km and p = Q.Qg/cm? (V ^ 2 x 10~^^eF) 
and best-fit KamLAND data P2| parameter values: Am^^; = 7.1 x 10~^eF^ 
and Osoi = 32.5°. 
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Figure 3: Neutrino oscillation probability Pe^e vs neutrino energy, when the 
neutrino flux pass through the center of the Earth with best-fit KamLAND 
data j22j parameter values: Aml^i = 7.1 x lO'^eV"^ and Ogoi = 32.5° 
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